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Abstract 

Transition rates for coherent muon-electron conversion in muonic atoms, (jlN — > 
eN, are computed for various types of muon number violating interactions. Atten- 
tion is paid to relativistic atomic effects, Coulomb distortion, finite nuclear size, and 
nucleon distributions. Discrepancies with previously published results are pointed 
out and explained. Results are presented for several elements of current and future 
experimental interest. 

1 Introduction 



In the standard model (SM) of electroweak interactions, leptons of different flavors do not 
mix because of vanishing neutrino masses. On the other hand, lepton flavor violation is 
predicted in various extensions of the SM. Hence, its observation would be direct evidence 
of physics beyond the SM. 

Because the muon is a relatively stable particle which can be abundantly produced, 
muon-number violating processes are of particular interest. Most recently, the follow- 
ing processes have been studied experimentally (for a review see e.g. [Jj, 3): // — > ej, 
muonium-antimuonium oscillations, and, the subject of the present work, muon-electron 
conversion in the field of a nucleus. Muon number violation has also been searched for in 
kaon, tau, and Z decays. 

■"■Talk given at the Workshop on Physics at the First Muon Collider and at the Front End of a Muon 
Collider, Fcrmilab, November 6-9, 1997. 
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A particularly sensitive search for muon number violation is the coherent conversion in a 
muonic atom, formed by muon and a nucleus. Here by coherent conversion we understand 
a process fi~N — > e~N in which the nucleus N remains in its initial state (up to recoil 
effects). The rate of the coherent conversion is enhanced with respect to processes with a 
nuclear excitation by a factor of the order of the number of nucleons. 

Since the original pioneering papers on the conversion theory || |], |5], there have been 
several theoretical efforts intended to provide an accurate description of coherent muon 
conversion. There are two kinds of theoretical issues: first, the short distance effects 
which are responsible for the muon number violation (presumably caused by some "new 
physics"); second, the long distance atomic physics of the muonic atom, which also de- 
termines the transition rate. 

The first group of problems has been studied in many extensions of the standard model 
(for a review see [H, |^, ||). In particular, in ref. || the rate of the coherent conversion 
fi~N — > e~N was calculated in a variety of gauge models. It was pointed out in that 
paper that in a large class of models the conversion can be much more probable than the 
decay /i — > erf. For example, there can be logarithmic enhancements of the form factors 
leading to conversion which are absent in the decay rate. Such logarithmic effects were 



also recently discussed in [10 



Weinberg and Feinberg || focused on an electromagnetic mechanism of transferring the 
energy yield to the nucleus. The structure of this electromagnetic interaction is richer 
than for the /i — > decay since the photon need not be on mass shell. Therefore, it 
is possible that conversion can occur even if /i — > &~f is forbidden for some reason. The 
matrix element for conversion contains monopole terms which do not contribute to the 
decay \i — > ej, because the longitudinal polarization states are possible only for virtual 
photons. In addition, there may be conversion amplitudes other than photon mediated 
processes. This makes the conversion on nuclei a particularly interesting process to study. 

The early theoretical studies of muon conversion into electrons on nuclei performed in 
H' H> HI are valid mainly for conversion on light nuclei. The developments before the 
year 1978 have been summarized in ||. In heavier atoms new effects become important: 
relativistic components of the muon wave function, Coulomb distortion of the outgoing 



electron, and the finite nuclear size. These were first addressed in 111. Nuclear effects 



were also analyzed, albeit in a non-relativistic approximation, in [12| and, more recently, 

mimm. 



Recently, we have undertaken a new calculation of the full relativistic atomic physics 
aspect of the conversion. In this talk our main results are summarized. The details of the 



calculation and a more extensive analysis will be given in a forthcoming paper [T5I . 



2 General description 

The muon-electron transition can occur via various mechanisms (see e.g. || [16], |l7l). To 
keep the description fairly general, it is convenient to write down a low energy effective 
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Hamiltonian for the fx — > e transition: 

T-l = -eOfi + h.c, (1) 

A a (q) (2) 



O = -V4 



not 



q 2 

7a {fEO ~ fMOlb) —o + icr af3— (/jWl + fElls) 

m 2 m 



G F 

+ -y|7a(« - &7s)^ Q , 

J a = u^u + Cdd^d. (3) 

In this equation m = is muon mass and a, b, /bo,i, /mo,i are dimensionless coupling 
constants. 

The part of the Hamiltonian containing the photon field A a describes the transition of 
a muon into an electron and an on-shell photon, fx — > ej, whose rate is T(fx — > e7) = 
otm(\f M i\ 2 + \ fEi\ 2 ) /2, which gives the branching ratio 5r(/i — > ej) = F(fx — > ej)/T(fx ez/gZ/^) 
96n 3 a (|/mi| 2 + l/sil 2 ) /(G^m 4 ). Due to gauge invariance, the transition vector and axial 
currents (terms proportional to Jeo, /mo) cannot contribute to the transition to an on-shell 
photon ||. We explicitly account for this by factoring out q 2 /m 2 . We further assume that 
all effective coupling constants introduced in the Hamiltonian are slowly varying functions 
of all external parameters (photon off-shellness, muon off-shellness etc.) and therefore 
can be considered as constants in the course of the calculation. 

To calculate the transition rate for the coherent conversion fxN — > eN we first average 
the effective Hamiltonian over the nucleus. The result of this average depends on muon 
and electron wave functions and the field of the nucleus: 

H int = (N\H\N), (N\N) = 1. (4) 

To calculate a matrix element of an arbitrary operator Q between two nuclei, we use the 
following approximation: 

(N\Q\N) = J d 3 r (^Zp p (r) (p\Q(r) \p) + (A - Z)p n (r) (n\Q(r) \n)j . 

In the above equation Z is the number of protons in the nucleus N and (A — Z) is 
the number of neutrons. Also, \p) and |n) denote states of a single proton and neutron 
respectively, with densities normalized as follows: 

d 3 rp p (r) = 1, J d 3 rp n (r) = 1. (5) 

To obtain matrix elements of the Hamiltonian, one should calculate the matrix element of 
the quark currents between two nucleons. A typical momentum transfered by the current 
is of the order of the muon mass, a small quantity compared to the mass of the nucleus. 
The time component of the soft current counts the number of constituent quarks in the 
nucleon. We therefore obtain 

(p\uj°u + c d dj°d\p) = 2 + c d , 

(n\u-f°u + c d d-f°d\n) = l + 2c d . (6) 
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The matrix element of spatial components of the current is proportional to the velocities 
of the constituents and is negligible in the present problem. 

Calculating the matrix element of the Hamiltonian with respect to the nucleus states, we 
arrive at an effective Hamiltonian for the coherent // — > e conversion in the field of the 
nucleus: 

Hmt = Hi + H 2 + H 3 , 
Hi = ej d 3 r$ e (r) 7o (f EQ - f M ol 5 ) ^ „{r)An Z p p {r) 

H 2 = ej d 3 rty e (r)ta af3 (f M1 + f Ell5 ) ^ „(r)F aP (r) 
H 3 = ^ J d 3 r* e (r) 7o (a - bj s )%(r) 

x (z(2 + c d ) Pp {r) + {A - Z)(l + 2c d ) Pn {r) ] j (7) 

In the above equations the fields \l/ e (r) and ^(r) stand for the second quantized operators, 
and the electromagnetic tensor F a p is understood as a classical electric field produced by 
the nucleus. The matrix element of ifi nt taken between appropriate initial and final states 
will give the amplitude for the coherent conversion. 

In our case the initial state is the muon in the IS* orbit around the nucleus; the final state 
is an electron with the energy equal to the energy of the initial muon. The corrections 
due to the nucleus recoil are small and we do not consider them here. 

In their pioneering work on muon-electron conversion Weinberg and Feinberg per- 
formed an approximate calculation of the coherent conversion rate. Using simple esti- 
mates, we would like to show when complete treatment of the problem requires going 
beyond the approximations used in Ref. 0. For this purpose we describe the scales 
relevant for the problem. The wave function of the muon bound in the lowest orbit is 
characterized by the Bohr radius a# = (aZm)" 1 . The radius of the nucleus Rn scales like 
m^R N ~ (Z/4) 1 / 3 . Evidently, the nucleus can be considered as point-like only if the Bohr 
radius is much larger than the radius of the nucleus a# 3> Rn- This implies Z <C 60. On 
the other hand, the relativistic corrections are governed by the parameter Za. Therefore, 
for high Z elements such as Pb, it is not clear a priori if the non-relativistic treatment 
of the muon bound state is sufficient. Moreover, another physical effect is governed by 
the same parameter Za. Consider the case when /i — > e conversion occurs due to an 
exchange of the photon with the nucleus. In this case the kinematics of the decay dictates 
that the virtuality of the photon is determined by the mass of the decaying muon. The 
process can be considered point-like, if this scale is much less than the Bohr radius. This 
implies that the photon-mediated conversion can be considered as a point-like process 
for Za 1. Therefore, for heavy nuclei such as Pb the consideration of the conversion 
process as point like is no longer valid. The appearance of the part H 2 in the effective 
Hamiltonian ifi nt reflects this observation: one notes that the /x — >• e transition current 
couples to the electric field of the nucleus, not to the proton density directly. 

Clearly, for light nuclei one can rely on the hierarchy of the scales and perform an approx- 
imate calculation of the rate. In this case the conversion rate will be proportional to the 
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square of the muon wave function at the origin and the square of the nucleus form-factor. 
Such an approximation was used in Ref. ||. When the charge of the nucleus grows, all 
scales relevant to the problem become comparable and the above approximations cannot 
be trusted. In this case a correct treatment of the problem requires solving the Dirac 
equation for both muon and electron wave functions in the field of the nucleus. 

To describe proton and neutron distributions in the nucleus we use two-parametric Fermi 
functions: 

Po 



P P (n)(r) 



I _|_ e ( r - r p(n))/«p(n) 



d rp p{n) (r) 



1. 



In the above equation r p(n ) ~ A 1 / 3 is the radius of the proton (neutron) fraction of the 
nucleus and a p i n ) is the thickness of the boundary of the proton (neutron) fraction. Precise 



values of these parameters depend on the nucleus and can be found in tables [18|, [19| . 
The transition rate for the coherent conversion is given by: 

u{pTN -> e~N) = uj conv = J2 \H[: t \ 2 . (9) 



Here the sum goes over all quantum numbers which the electron in the final state can 
have in addition to energy. The muon bound state wave function is normalized to unity. 
The electron wave function of the continuous spectrum is normalized such that: 



d 3 n/^, x/ (r)^E,x(r) = 5 x , x ,2n5(E' - E) 



(10) 



where A denotes a set of all discrete quantum numbers which electron obeys in addition 
to energy. 



It is convenient to use the expression for the wave functions as suggested in [2(| : 



1> 



g [) ( r )xT 

j/(*)( r ) x ™ 



In this equation m is the eigenvalue of the operator J z , J = L + S; —k is the eigenvalue of 
the operator L ■ a + 1 (see e.g. |2(|); xT is an orthonormalized spinor. For the muon wave 
function we use the lowest energy bound state wave function; this suggests that j = 1/2, 
m = ±1/2, k = —1. We will distinguish muon and electron radial wave functions by a 
subscript. 



A calculation of the explicit expression for the transition rate uj c 
We find 



is now straightforward. 



0J r 



drr 



drr 



fE (r)(g e g^+f e f^) + ^(g e 4 + f e <7 M )^ 

f m dV 

fM0{r){g e g fi +f e f^) + —(9e f/j, + fe9^)~^: 



ho{r) 
fMo(r) 



-f E0 ^^P P (r) + [Z{2 + c d )p p (r) + N(l + 2c d )p n (r)} 

-fuo^^Ppir) + -|6 [Z{2 + c d ) Pp {r) + N{1 + 2c d )p n {r)} 
m A a/2 



(12) 
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where V(r) = — eA°(r) is the muon potential energy in the field of the nucleus. 

The above equation provides a general expression for the rate of the reaction p~N —>■ e~N 
which we will use for numerical analysis in the next section. The radial wave functions 
/~ and are obtained by solving the Dirac equation. Inside the nucleus and close to it 
these solution are found numerically. At large distances we match the numerical solutions 
to the exact Coulomb wave functions. 



3 Numerical analysis 



In general, the following integrals are needed for the description of the transition rate: 



n 
/ 3 



4nZa 

m 2 
AnZa 



drr 2 p p (r)g g 



fi tie i 



drr 2 p n (r)g g, 



tie i 



1 

rn 



,dV(r) 



dr 



(j, J e i 



1 2 



l 2 ~ 
1 

m 



AnZa 

m 2 
AnZa 



drr*p p (r)f f e 



drr 2 p n (r)f^ f e 



drr — / g t 



dr 



H tie 



(13) 



Tables with the values of these integrals for various elements will be given in [15]. Using 
those values, the conversion rate is 



° 1U V^conv i 10 com 



sec 



(14) 



(i) 



f E0 I p - -^-^— a (z(2 + c d )I p + N(l + 2c d )I n ) + f m I^ 



(2) 



y/2AnZa 
G F m 2 



v^AnZa K - 



b(z(2 + c d )I p + N(l + 2c d )I n ) + f E1 I u 



where 



/ p = -(/f+/ 2 p ), /„ = 



'34 



h + h- 



(15) 



(16) 



All dimensional parameters in the above equations are expressed in fermi. For proton and 
neutron distributions we use the two-parameter Fermi distribution (§) with parameters 
taken from Ref. 

For the present application we neglect the effects of the vacuum polarization; such ap- 
proximation is justified in view of much bigger errors in wave function integrals induced 
directly by the nuclear distribution uncertainties. 



4 Branching ratio for fi-e conversion in a specific model 

We consider here one specific model which predicts lepton flavor violation, a supersym- 
metric grand-unified theory discussed in Ref. |T7| . 
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The analysis of Ref. |17|] implies that magnetic couplings /mi and f E i (cf. Eq.(2)) 



are significantly enhanced in a large region of the parameter space of those models in 
comparison with other couplings in the effective Hamiltonian. Therefore, it is reasonable 
to neglect all other couplings in the effective Hamiltonian and to analyze the dependence 
of uj conv /u capt on the choice of target. 

With these approximations the ratio of the conversion rate to the capture rate becomes 



^capt 



3-10 12 (\f E1 \ 2 + \f M1 \ 2 )B(A,Z), (17) 



I 2 

with B(A,Z) = 10 n f -. (18) 

Ucaptl sec- 1 

We compare the above formula with the Weinberg-Feinberg approximation, in which the 
rate is obtained by replacing B(A, Z) in ([17|) by B WF (A, Z): 

B WF (A, Z) = 8a 5 mZ± s ZF* 1 1 (19) 

^capt O • 1U 

In this formula Z q q denotes the "effective Z" obtained by averaging the muon wave 
function over the nuclear density (see e.g. ||12||), and F p is the formfactor describing the 
charge distribution, given by 

f d 3 ro(r) sinmr 
p ~ Jd 3 rp(r) ■ [ > 

In table |l] we show the B and Bwf coefficients for three elements. We conclude that 
in this particular model the ratio of the coherent conversion rate to the capture rate, as 
described by our B(A,Z), does not change significantly with changing the target. The 
WF approximation tends to slightly overestimate the conversion rate. In ref. flT 



various 



corrections to the WF approximations have been studied; in that approach B WF should 
be replaced by B s = CiC 2 C 3 B WF , where the correction factors Ci^z are listed in table III 
in ||11||. We agree with these results for non-photonic mechanisms of conversion. However, 



for the photonic case we are considering here, our results differ from [Tl[] even stronger 
than from the WF approximation, as can be seen in table [I]. This is because the effect 
of the non-locality in the interaction with the electric field was not taken into account in 



1 1 1 . It is especially important for heavy elements. 



It is instructive to compare the conversion rate to the branching ratio for fi — > in this 
model: 

Briii -> ej) 96vr 3 a 1 428 

(21) 



,Ma P t G\m^ 3 • 10 12 B(A, Z) B{A, Z) ' 



This ratio varies from 389 for 27 A1 to 238 for 48 Ti, and increases again for heavy elements 
to 342 for 208 Pb. 
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Table 1: Comparison of our results for the coefficient B(A,Z) in the rate formula ( |i~7|) 
with the Weinberg- Feinberg approximation Bwp(A, Z) and with Bg{A, Z) obtained from 
ref. |TT[. Z q q is taken from |L^] and the capture rates u) capt from [pTJ] ; F p is computed 
using (H). 



Element 


B(A,Z) 


B WF (A, Z) 




Z eff 


F P 


Ucapt, [10 6 /sec 




Al 


1.1(1) 


1.2 


1.3 


11.62 


0.63 


0.7 


Ti 


1.8 


2.0 


2.2 


17.61 


0.53 


2.6 


Pb 


1.25(15) 


1.6 


2.2 


33.81 


0.15 


13.0 



5 Bounds on the muon number violating couplings 
from SINDRUM 

Using updated results for the atomic part of the theoretical description of the coherent 
muon electron conversion, we reconsider the analysis performed in Ref. p2 1 . In that paper 



the upper limit for the ratio tOconv/^capt has been reported based on the measurement using 
208 Pb. These results were combined with previous measurement on 48 Ti target to obtain 
the bounds on muon-number violating coupling constants. In that analysis it was assumed 
that the effective Hamiltonian does not contain a piece which corresponds to the photon 
mediated conversion. 



As a first step we have to switch to the model and notations of Ref. ||22|| . The model, 
considered in Ref. [P2"| , corresponds to the exchange of a heavy vector or scalar boson 
which mediates \i — > e transition. For this reason, we should equate /#o,i and /mo,i to zero 
in our general expression for conversion rate, eq. (|12|)- Also, to switch to their notations, 
one should substitute b = a, aca — > (g v — g v )/2, and a — > (g v + g v )/2 in that equation. 
Constants gy,gy then parameterize the coupling of the heavy gauge boson to isoscalar 
and isovector parts of the quark current respectively. 



Then, one gets for the conversion rate: 



A av \ 3A 




(22) 



In this equation we use the notation (cf. Eq. (|T^)): £ = I p /I n - Using the bounds on the 
\l — > e branching ratios measured with Ti and Pb targets, we derive the new bounds for 
the coupling constants g v , g v : 

\g v \ < 8±2- 10" 7 , |^| <40±13-10~ 6 . (23) 

The variation in the boundaries shown above is the estimate of the uncertainty of the 
result induced by uncertainties in the input nuclear parameters. 

These new bounds should be compared with the values 

|<7°| <3.9-10- 7 , |^| O.7-10- 6 . (24) 



S 



which were given in Ref. |22| using the results of Ref. |TT . 

There are two reasons for such large differences. As we have mentioned already, in Ref. 
P2] the results of Ref. [O] have been used to interpret experimental data. We note in 



this respect, that Table I of ref. |TT| does not provide correct results - the results for 
Uconv/u capture quoted there are about a factor of 2 too large. Second, in the case of a 
Pb target the difference in proton and nucleon distributions becomes quite noticeable. 



However, the numbers in Table I of [TTJ were obtained using identical proton and neutron 



distributions. These two reasons conspire to give quite a large discrepancy. 

We note, however, that the values for the bounds quoted above (especially for g\?) are 
very sensitive to the input parameters. This happens, because the bounds are derived 
from a system of linear equations in which the coefficients of g\ there are much smaller 
than all other parameters. Therefore even a small variation in the input parameters can 
generate quite substantial change in the resulting value for g\. 
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